Schwinger Pair Production at Finite Temperature in Scalar QED 

Sang Pyo KinE 

Department of Physics, Kunsan National University, Kunsan 573-701, Korea and 
Asia Pacific Center for Theoretical Physics, Pohang 790-784, Korea 

Hyun Kyu Le^ 

Department of Physics and BK21 Division of Advanced Research 
and Education in Physics, Hanyang University, Seoul 133-792 and 
Asia Pacific Center for Theoretical Physics, Pohang 790-784, Korea 

(Dated:) 

In scalar QED we study the Schwinger pair production from an initial ensemble of charged bosons 
when an electric field is turned on for a finite period together with or without a constant magnetic 
field. The scalar QED Hamiltonian depends on time through the electric field, which causes the 
initial ensemble of bosons to evolve out of equilibrium. Using the Liouville-von Neumann method 
for the density operator and quantum states for each momentum mode, we calculate the Schwinger 
pair-production rate at finite temperature, which is the pair-production rate from the vacuum times 
a thermal factor of the Bose-Einstein distribution. 

PACS numbers: 12.20.-m, 13.40.-f, ll.lO.Wx 



I. INTRODUCTION 

Vacuum polarization and pair production has been known more than seven decades ago Schwinger used the 
proper-time method to find the exact one-loop effective action for a uniform electromagnetic field Since then 
vacuum polarization and Schwinger pair production has been investigated more or less from a theoretical view point 
(for a review on recent development of strong QED, see Ref. Q). However, recently the study of strong QED is more 
physically motivated than ever since it will have many interesting physical applications in the near future. In the near 
future terrestrial experiments, the x-ray free electron lasers from Linac Coherent Light Source at SLAC and the TeV 
Energy Superconducting Linear Accelerator at DESY are expected to achieve electromagnetic fields near the QED 
regime 0, H, Q. In astrophysics, most neutron stars are expected to have magnetic fields comparable to the critical 
strength Be = m?c^ /eh (4.4 x lO^^G for electrons) and furthermore some magnetars are known to have magnetic fields 
a few order above the critical strength. It is interesting that the ratio of the pair-production rate in the presence of 
parallel electric and magnetic fields to that in the pure electric field is Af{E, B)/J\f{E, 0) — {irqE /qE) x coth{7TqB / qE) . 
This implies that the pair-production rate is enhanced by the strong magnetic field, though the leading behavior is 
still determined by the electric field as /\f{E,0) — {iqE)^l'^e~^™ /'?-^/(47r'^m). In these compact stars pairs would be 
produced at non-zero temperature. 

Another interesting issue from the view point of strong QED and astrophysics is strange quark stars, hypothetical 
astrophysical objects. Witten proposed that strange quark matter niight be a stable state of nuclear matter at 
ultrahigh density, which was numerically studied by Farhi and Jaffe Based on this hypothesis strange quark 
stars were proposed as an alternative to neutron stars 0] . The charge neutrality and chemical equilibrium of strange 
quark stars would require a net amount of electrons, which are bounded by the Coulomb attraction from the positive 
core and form an electrosphere on the surface. The Coulomb barrier from the electrosphere can generate an electric 
field as large as 5 x which is stronger than the critical strength Ec = rri^c^Jeh (1.3 x IQ^^V/cm). The 

Coulomb barrier may be a source of electron and positron pairs at finite temperature . Therefore the Schwinger 
pair production at finite temperature would be interesting for such astrophysical objects. 

It was Dittrich who first found the QED effective action at finite temperature in the presence of a constant magnetic 
field [lOj. Since then this has motivated the study of QED at finite temperature under various configurations of 
electromagnetic fields and thereby the Schwinger pair production at finite temperature, which has been an issue of 
debate depending on formalisms employed. The QED effective action for a constant magnetic field was obtained 
at finite temperature and density [11]. Loewe and Rojas used the Schwinger proper-time method to obtain the 
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effective action in the presence of both constant magnetic and electric fields and found an imaginary part that has 
contributions from vacuum as well as thermal fluctuations, leading to a thermal enhancement of pair production 
at high temperature [l^. On the other hand, Elmfors and Skagerstam used the real-time formalism to obtain the 
QED effective action in the presence of both magnetic and electric fields, in which there is no imaginary part and 
thereby no pair production [l^]. Contrary to this, Ganguly et al found in the finite temperature field theory that 
the pair-production rate increases at high temperature [ij]. Also in the functional Schrodinger picture Hallin and 
Liljenberg found the density matrix for QED when the electric field is turned on only for a finite period and thereby 
calculated the pair production at finite temperature when the system evolves from an initial ensemble ,15j . Recently 
Gies used the imaginary-time formalism to obtain the QED effective action in the presence of both constant electric 
and magnetic fields, which does not have any imaginary part at one-loop but nonzero imaginary part at two- loop p^ . 

The main purpose of this paper is to study the Schwinger pair production from an initial ensemble of charged bosons 
in scalar QED when an electric field is turned on only for a finite period together with or without a constant magnetic 
field. A Sauter-type electric field E{t) = Eq sech^ (t/r) provides an analytical model p7|], which acts effectively for 
— T < t < T. Thepair production from vacuum by this profile of electric field extended in a finite time or space 
has been studied [11, E, H [HI, [H [H 111 . Note that as the gauge potential = (0, 0, 0, -£;or tanh(t/r)) in the 
time-dependent gauge depends on time explicitly, so does the Hamiltonian. As the evolution operator e"''-^ in the 
Euclidean time is not the density operator, the imaginary-time method based on the Matsubara frequency may not 
be used directly for general time-dependent cases. In fact, the initial ensemble evolves out of equilibrium due to the 
time-dependent Hamiltonian, which requires a nonequilibrium quantum field theory. 

In this paper we shall use a nonequilibrium quantum field theory based on the Liouville-von Neumann equation, 
which was o rig inally introduced to find the density operator but led to exact quantum states of time-dependent 
systems [25l . [26| . The scalar QED Hamiltonian is quadratic in the field and momentum and each momentum mode 
has a time-dependent frequency. The time-dependent creation and annihilation operators that satisfy the Liouville- 
von Neumann equation provide not only the density operator but also exact quantum states. We find that the number 
of pairs produced by the electric field is proportional to the initial Bose-Einstein distribution. From this scalar QED 
model we may conclude that the Schwinger pair-production rate by a time-dependent electric field is enhanced by a 
thermal factor of the initial Bose-Einstein distribution. The Liouville-von Neumann method can be compared with 
other methods used to handle the scalar QED Hamiltonian. Hallin and Liljenberg used the Schrodinger picture to 
find the wave functional of the Schrodinger equation and the density matrix lE\ , whereas the authors of Refs. [13, [ll] 
used the Heisenberg picture to find the time evolution of the creation and annihilation operators for each momentum 
of the field. It should be remarked that the time-dependent creation and annihilation operators of the Liouville-von 
Neumann equation provide not only the density operator but also exact quantum states for each momentum. 

The organization of this paper is as follows. In Sec. II, we explain the canonical quantization method based on 
the Liouville-von Neumann equation and apply it to scalar QED in the presence of an electric field. In Sec. Ill, we 
calculate the Schwinger pair-production rate at finite temperature in the presence of Sauter-type electric field together 
with or without a constant magnetic field. 



II. CANONICAL QUANTIZATION OF SCALAR QED 



Scalar QED for a charged boson with mass m and charge q under an external electromagnetic field with the gauge 
field Af^ is described by the Lagrangian [in units with h = c = I and with metric signature (+, — , — , — )] 



L = I (tx 
where the covariant derivative is given by 



(1) 



(2) 



For instance, the electric field in the z-direction may be described by the time-dependent gauge potential = 
(0,0,0, ^z). Then the Klein-Gordon equation from the Lagrangian ([T]) takes the form 

~dl- id, + iqA,it)f + m2]$ = 0. (3) 

The corresponding Hamiltonian is given by 



H = J (fix U*n + <^>* (^-dl - {d, + iqA^f + m^) $ 



(4) 
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where 11 = $* and 11* = $ are the conjugate momenta. 

Following Ref . ^] , the fields and momenta are decomposed into Fourier modes as 

n(x,i) = fm^le-''^-^ = /[rfk]^k 



(5) 



where [dk] = fi^k/(27r)^ and 



n*(x,i) ^ / [dk]0ke"" 



ik-x 



The Hamiltonian is then given by 



where 



H= [dk] 



^kW = (^'^ + qMt)f + ki + m^. 

Here, 4'\5.i4't. ^^'^ "'kj^'k ^'''^ Schrodinger operators. Each Fourier mode obeys the equation 



(6) 



(7) 



(8) 



(9) 



Upon canonical quantization, we impose the equal-time commutation relation [$(x', f), n(x, t)] — iS{x' — x) and 
[$*(x', t), n*(x, i)] = iS{x' — x), from which follow the commutation relations 



/•k^^k] - ^(27r)35(k' - k), 
ru',K] = *(27r)3<5(k' - k), 



(10) 



but all the other commutators vanish. We may introduce the time-dependent annihilation and creation operators to 
quantize the fields as 



and similarly the momenta as 



<i>(x,t) 
<i>*(x,t) 

n(x,t) 
n*(x,t) 



[dk] 
[dk] 

[dk] 

[dk] 



^kW^kW + ^kW^kW 

^kW&lW+^k(i)akW 



-zk. x 



-ikx 



^■ik. x 



(11) 



(12) 



A few comments are in order. First, note that and aj^. correspond to the particle and 6k and to the antiparticle 
and that $ annihilates particles but creates antiparticles. Second, it can be shown that the field quantization includes 
static fields as a specific case by choosing (pk = e~'"''*/\/2'jJk- Finally, by imposing the Wronskian conditions 



we have the standard equal-time commutation relations 

[ak'(i),4W] = S{k' - k), [6k'(i),4W] = - k). 
All the other commutators vanish. 



(13) 



(14) 
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By inverting the Fourier modes in Eqs. (jlip and (|12p . we find two independent sets of time-dependent annihilation 
and creation operators [261] 



and 



(15) 



(16) 



It is remarkable that when iy9k satisfies the mode equation ^ , the time-dependent operators ([T5t and (|16p satisfy the 
Liouville-von Neumann equation, for instance, 



.5ak(t) 

* dt 



M),H{t)]^0. 



This implies that we may define the density operator for particles 



/ \ 1 r / t / \ / \ 1 

PaM = ^*^^Pr fc^V"kWak(i) + 2 



and for antiparticles 



(17) 



(18) 



(19) 



Here, Zk = Tr(e-"''(''k(*)«kW-Hi/2)/fcT) ^ rp^(g-<^k(f>i,(t)6>c(t)-Hi/2)/fcr-)^ -y^^ ^^^^^ ^^^^^ ^ ^j^^ temperature of an initial 
ensemble of the heat bath of particles before the onset of electric fields and that the number of produced pairs should 
be small not to change the temperature of heat bath. 

We now turn to pair production by electric fields. The mode-decomposed Hamiltonian ([7]) is a collection of time- 
dependent harmonic oscillators. The equation of motion in the language of quantum mechanics, becomes the 
super-barrier transmission over a potential barrier in the time dimension. Pair production may be explained in such 
a way that the positive-frequency solution at the past infinity scatters not only to the positive-frequency solution but 
also to the negative-frequency solution at the future infinity. Then the magnitude square of the ratio of the coefficient 
of the negative-frequency solution at the future infinity to the coefficient of the positive-frequency solution at the 
past infinity is the pair-production rate [2^. [soj. However, we shall employ the Liouville-von Neumann method, in 
which we first find the operators known as the Lewis-Riesenfeld or invariant operators that satisfy the Liouville-von 
Neumann equation [2^1 and then construct not only the exact quantum states but also the density operator. The 
time-dependent vacuum is defined as 



-u;i,{bl{t)b^,{t) + l/2)/kT\ 



ak(0|0;i) =6k(t)|0;<) =0, (Vk) 
and the multi-particle and antiparticle states as 



4rw cw 



t"2. 



Pki • • • ;nk2 • 



■,t) = 



/ni\ 



'n2\ 



\Q-t). 



The time-dependent vacuum ()20|1 is the zero-particle and antiparticle state of the number operators 

Na^{t) = al{t)ak{t), N,^it) = &],(i)6k(0- 



(20) 



(21) 



(22) 



It should be remarked that the Fock states are exact quantum states of the time-dependent Hamiltonian up to some 
time-dependent phase factors a la the Liouville-von Neumann method [l^, [2^ . 

To compare with the standard quantization of field, we consider the limiting case of no electromagnetic field. Then 
our definition of the time-dependent annihilation operator differs from the standard one by an overall time-dependent 
phase factor: 



ak(t) 



(23) 
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for the choice of (pk = e~'"''*/-\/2tJk. The phase factor, which is necessary for the Liouville-von Neumann equation, 
does not change the number operators -/Va^(t) = Na^. and -/V&k(i) = Nb^. In this sense the Liouville-von Neumann 
method includes the standard quantization as a specific case. The method turns out to be convenient and powerful 
for time-dependent quantum fields such as scalar QED in the presence of time-dependent electromagnetic fields. The 
number of pairs with a certain momentum k produced by an external field is the number of particles or antiparticles 
of the out-vacuum at t = oo contained in the in-vacuum at t = —oo pol. [30l| 

nk = (0;in|iV,Ji = oo)|0;in) = (0; in|7Vb^(t = cx))|0; t = in). (24) 

III. PAIR PRODUCTION AT FINITE TEMPERATURE 

We first study the Schwinger pair production by a uniform field E{t) — Eosecil^ (t / t) without a magnetic field at 
finite temperature and then study the case with a constant magnetic field B. When E • B 7^ we can find a Lorentz 
frame where the electric and magnetic fields are parallel to each other, for instance, in the z-direction. We choose 
is a Sauter-type gauge potential = {0,—By/2,Bx/2,—EoTtanh{t/T)). It is shown that pair production by the 
Sauter-type potential is characterized by two parameters [22I ] 

et = Ot = J. (25) 

The pair production is energetically favored for < 1 and the adiabaticity parameter 5t determines the pair-production 
rate. 

A. Pure Electric Field 

The mode equation for the auxiliary field variable takes the form 

(p^,{t) + ujl{t)^^,{t) = (26) 

where 

U}l{t) = [k^ - qEQTtimh{-)^ +'k\+ (27) 
Each mode has two asymptotic frequencies ai t — =pcx3 

= y^(fc.-gi?oT)2+ki+m2. (28) 

The solution with the asymptotic form 

(plf (t) = ^= (29) 



24-) 



at t = —00 is given by [H, [T^ 



^k(i) = , ^ (1 - zY'^+^^z-"^-''^F{a,,, /3k; 7k; z), (30) 



24 ^e'^^k ' 

where i<"(ak, /3k; 7k; -s^) is the hypergeometric function and 

z = -2e2*/- 



A = ^j{qE,T^Y-\ (31) 
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and 



1 i 



1 i . (_) 

k 

(-) 



1 i . 



2 2^~'' 

7k = 1 - 



(+) 



-2A), 
2A), 



(32) 



Note that the solution l[30|) satisfies the Wronskian condition (fT3|) for any time. At the other asymptotic hmit t = oo, 
the solution becomes 



where 



„(+), 



2w 



(+) 



(33) 



(34) 



and 



,(+) 



r(7k)r(/3k - Qk) 
^4-^r(/3k)r(7k-«k)' 



,,(+) 



,(+) 



r(7k)r(ak - /3k) 

^^wr(ak)r(7k-/3k)- 



(35) 



As the asymptotic form (|33p is approximately valid for t > t, the analysis below is a good approximation after r and 
exact at t = oo. 

The operators in Eqs. (fT5|) and (fTB|) obtained by substituting the asymptotic solution (I29p define the time-dependent 
annihilation and creation operators flk'i'^k" ^'^^ ^Ifi^k" ^^"^ in-vacuum while those obtained by the asymptotic 
solution define the time-dependent annihilation and creation operators a^^*,aj^°"' and 6^^*, 6]^°"* for the out- 
vacuum. The mode solution starting fi'om the initial asymptotic solution (|29p evolves to another asymptotic form 
([55)1 . This leads to a Bogoliubov transformation between the in-vacuum operators and the out- vacuum operators 



^out 



A in I D* i.tiii 



^k6i:' 



whose coefficients satisfy the relation 



(36) 



(37) 



The in-vacuum is the superposition of the out-vacuum particles. The number operator of the in-vacuum is iV^" = 



.tin. 



for particles and N^^ — h^^U^ for the antiparticlcs. Similarly, the particle number operator of the out-vacuum 

Hence the number of created pairs by 



is Ar°"' = a!°"*"°"* 



■ Qk "'k "-k 

the electric field from the vacuum is 



and the antiparticle number operator is N'^^ 



out _ ^fout^out 



where 



nk = (0;in|A^-"|0;in) = (0; in|iV,™*|0; in) = \B^.\\ 

cosh[7r(2A Tw^+^ - tJ~^)/2] cosh[7r(2A + tJ~'^ - ^['^^2] 
cosh[7r(2A + tJ^^ + Tuj^^^)/2] cosh[7r(2A - tJ^ ^ - tJ^^)/2] ' 

The number of pairs (|39p approximately is 



|SkP«e- 



e k 



(38) 



(39) 



(40) 
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where 5^°'' is the leading-order contribution of the WKB instanton action in scalar QED [22 
The density operator takes the asymptotic form for particles 



(-) 



and the density operator (fTO]) for antiparticles 



1 



■ exp 



kT 



Hence the initial ensemble of particles has the Bose-Einstcin distribution 

1 



/k(T)^Tr(pi^W^;) 



e4 VfcT_ 



Finally, the number of pairs produced by the electric field from the initial thermal ensemble is 



nk(i?,r) = Tr(pi,"iV-*) - /k = \B^{E)\H2UT) + 1). 



(41) 



(42) 



(43) 



(44) 



The thermal factor can be written as 2/k + 1 = coth(a;^ '/2fcT). We have thus found that the Schwinger pair- 
production rate at finite temperature is enhanced by the thermal factor /k. In the zero-temperature limit we have 
/k = and recover the Schwinger pair-production rate from the vacuum. The Schwinger pair-production rate per 
unit volume and per unit time is given by summing over all momenta 



Af{E,T) = j[dk\n^{E,T). (45) 

As the solution ([30]) approximately approaches the asymptotic form (I29p for t « — r and (j33|) for t Ki t, the pair- 
production rates (j44p and (I45p are approximately valid after t = t. 

In the asymptotic region where the electric field is turned off, the out-vacuum is well-defined and the excited states 
correspond to a number of particles and/or antiparticles. In this region the pair-production rate (HH) or ([151) 
the interpretation of created pairs. Now a question may be raised how to interpret Eqs. (|44p or (|45|) while the electric 
field is acting. In terms of the in- vacuum operators we may express the time-dependent annihilation operator at any 
time as 



ak(0 = Au{t)4 



fok(i) = Au{t)b^ + B^^it} 



k ' 
'k ' 



where 



(46) 



(47) 



Here (ysj^ is the solution at t = — oo. Then we can obtain 



nu{E,T;t) ^ Tr pj^W^Ji) - /k - \BUE;t)\^i2MT) + 1) 



(48) 

It is tempting to interpret (|48p as the number of particles at time t. However, the concept of particle can have a 
meaning only in the adiabatic limit where 



(49) 



V2wk(i) 

is an approximate solution to the mode equation (^5]) . The solution is approximately good for a large momentum 
since Wk 3> |wk|, pk|- For a small momentum the change of the gauge potential and the electric field should be small 
to guarantee the validity of the adiabatic solution. In the adiabatic case the solution with the initial asymptotic form 
evolves approximately to the form 



^k(i)« A(tX(i) + Bk(iX*W 



(50) 



In the other case where the adiabatic approximation breaks down, we cannot have a concept of particle for such 
momentum modes but Eq. (|48p may be interpreted as a distribution. 
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B. Electric and Magnetic Fields 



We now consider the case of a uniform field E{t) — i?osech^(t/T) together with a parallel magnetic field B. 
gauge potential — {0, —By /2, Bx/2, ~Eqt ta,nh{t/T)). The mode equation takes the form 



where 



<iPk{t,x,y) +ujl{t)<fik{t,x,y) = 0, 
LuKt) = H±+ [k^ - qi;oTtanh(--)j + 



The 
(51) 

(52) 



Here H± describes the two-dimensional motion transverse to the magnetic field 



^qB \2 



Hi_ = -[dx-i—y] -[dy + i—x 



qB 



It has the eigenvalue 

H±ipn{x.,y) ^ qB{2n + l)ipn{x,y), 
whose orthonormal eigenfunctions were given in Rcf. [31]. Then the mode (pk{t,x,y) — ipnk{t)(pn{x,y) satisfies 

>fnk{t) + Ujlk{t)fnk{t) = 0, 

where 



2 



[k^ - qEoTta,iih{-)j +qB{2n+l 



(53) 

(54) 
(55) 

(56) 



The two asymptotic frequencies at t = ^oo are given by 

i^il' = Vikz + qEorY + qB{2n + 1) + m2, 

^nk 



<r = ^{k,-qEoTY+qB{2n + l)+^^ 



The mode has the solution of the same form as Eq. ([50)1 

, , 1 

<-Pnk\t) = 

with parameters replaced by 



(1 _ zY/^+'^Z~"'^'^>^''^F{a^k.Pnk-nnk;z), 



(57) 



(58) 





1 


* ( (-) 




^nk 


2 


(ink = 


1 
2 


« 1 (-) 




Ink = 


1 







The Bogoliubov transformations for particles and antiparticles are given by 

„out A in I p* Lfin 

;,out _ /I _ J,in ID* ^tin 



LOUt A tm I r>* „Ti" 



(59) 



(60) 



where 



B„ 



^ifc^ r(7„fc)r(/3nfc - o;„fc) 



'^iV r(7»fc)r(anfc - /3„fc) 



(61) 
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where 

/».(T)=Tr(pir,iV^„,)^ ^^,_,;^^_^ . (62) 

Finally, as in the pure electric field case, we find the number of pairs produced from the initial thermal ensemble 

n„fc(£;,B,T) = Tr(pj;>,7V™*) - U = \B,,k{E, B)\^(2U{T) + 1). (63) 

In the case of a pure magnetic field, ank = and thereby Bnk = as r(a„fc = 0) = oo and Ank — e^*'^""'' becomes a 
pure phase factor. This confirms the fact that pure magnetic fields cannot produce pairs of particles and antiparticles 
even at finite temperature. The total number of pairs per unit volume and per unit time is given by summing over 
the Landau levels and the longitudinal momentum 

Af{E,B,T)^-^Y. j dK\Bnk{E,B)\^{2fnk{T) + l). (64) 
Here [qB) /{2ti) is the number of Landau levels and another factor \/{2ti) is from the fc^ integration. 

IV. CONCLUSION 

In this paper we have studied the pair production of charged bosons at finite temperature by an electric field 
together with or without a constant magnetic field. As an analytically solvable model, the electric field of Sauter- 
type is used, which is effectively turned on for a finite period of time. Upon mode-decomposition the scalar QED 
Hamiltonian becomes a collection of time-dependent oscillators. We used the Liouville-von Neumann method to find 
the exact quantum states as well as the density operator for each momentum. The time-dependent gauge potential 
for the Sauter-type electric field has two asymptotic regions where the in-vacuum and out-vacuum are well-defined. 
The in-vacuum, which is annulled by the particle and antiparticle annihilation operators for each momentum before 
the onset of the electric field, contains a number of particles and antiparticles of the same momentum after the electric 
field is turned on and thus leads to pair production by the electric field. We have found that the number of created 
pairs for each momentum at the finite temperature is nnk{E, B ,T) = ri„fc(i?, i?, O)coth(a;^ '/2fcT), the number of 
pairs produced from the vacuum times a thermal factor from the initial Bose-Einstein distribution. In the limit of 
zero temperature the number of created pairs reduces to that from the in-vacuum as expected. We also confirm that 
the pure magnetic field does not produce pairs even at finite temperature at one-loop level. We may thus conclude that 
the Schwinger pair-production rate is indeed enhanced by the thermal effect given by the Bose-Einstein distribution. 

An interesting problem not handled in this paper is pair production of fermions at finite temperature by time- 
dependent electric fields. The fermion pair production by an arbitrary time-dependent electric field was studied in 
Refs. [S^ll^. Also the fermion pair-production rate by the Sauter-type electric field was found in the WKB instanton 
action method (2^ and in the worldlinc instanton method [23l [2^. However, to extend fermion pair production to 
the finite temperature case, we need a density operator for the time-dependent fermion system. For instance, the 
time-dependent annihilation and creation operators in Ref. [s^ l that satisfy the Liouville-von Neumann equation may 
be used to calculate the fermion pair-production rate at finite temperature. Another interesting problem is the back 
reaction of the produced pairs to the initial electromagnetic field. The back reaction problem in scalar and spinor 
QED was studied in Ref. [11]. As produced pairs can generate an additional electromagnetic field to the external 
field, the produced pairs provide a source term to the Klein-Gordon equation or the Dirac equation. The back reaction 
problem at finite temperature in scalar and spinor QED will be addressed in a future publication. 
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